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These lectures start with an elenientary introduction to the subject of 
magnetic monopoles which should be accesible from any physics back- 
ground. In the Weinberg-Salam model of electroweak interactions, mag- 
netic monopoles appear at the ends of a type of non-topological vortices 
called electroweak strings. These will also be discussed, as well as recent 
simulations of their formation during a phase transition which indicate 
that, in the (unphysical) range of parameters in which the strings are clas- 
sically stable, they can form with a density comparable to topological vor- 
tices. 



1. Introduction 

Last year marked the 50th anniversary of one of P.A.M. Dirac's most pro- 
found and famous papers - on magnetic monopoles. While his 1931 paper 
is usually considered the official birthday of magnetic monopoles, his 1948 
paper is where he really showed that it was possible to have a consistent 
quantum theory of magnetic poles in conjunction with electric charges, and 
described the interaction between them [1]. 

A standard reference on magnetic monopoles is Preskill's lectures in 
the 1985 Les Houches school [2]. Here I have tried to present magnetic 
monopoles and vortices in a way that makes them accessible to physicists 
who are not so familiar with the language of high energy physics, in par- 
ticular with e.g. non-abelian gauge theories. As a result, these lectures are 
much less technical. 
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The theme of this school is the use of topological defects as a tool 
to understand the dynamics of phase transitions out of equilibrium. It 
turns out that, since magnetic monopolcs and clcctroweak vortices are non- 
topological in the Weinberg-Salam model, their study can be particularly 
interesting in order to understand the role of gauge fields in defect formation 
during phase transitions. 

2. The elusive monopole 

The problem with magnetic monopolcs is well known in the context of 
Maxwell's equations. If the electromagnetic field is described by the vector 
potential A, then 

V-B = (1) 

so there can be no sources or sinks for the magnetic field. 

It would be fair to say that the experimental evidence for the existence 
of magnetic monopoles is not good. The 1998 Review of Particle Properties 
by the Particle Data Group [3] shows the result of monopole searches in 
particle accelerators: not a monopole in sight. On the other hand, cosmic 
ray searches have essentially only one event for which there seems to be 
no alternative explanation, observed by Cabrera in Stanford in 1982 [4]. 
The Cabrera detector, like many others that failed to find anything before 
or after it, consisted of a superconducting ring where a persistent current 
was monitored for a long time; in [4] the loop had an area of 20 cm^ and 
was monitored for a total of 151 days. During this time a single event was 
recorded which could be ascribed to a magnetically charged particle with 
one Dirac unit of 27rhc/e. A magnetic monopole. 

Under the circumstances, the experimental evidence is neatly summa- 
rized by the sentence (whose author is unfortunately unknown to me) "It 
is not clear that nobody has ever seen a magnetic monopole; what is clear 
is that nobody has ever seen two". And yet since Dirac's seminal work there 
have been over three thousand papers in the literature about magnetic 
monopoles! [5] 

The reasons behind this fascination with monopoles have evolved with 
time, but they are basically three: 

- the existence of monopoles would explain the quantisation of electric 
charge (for which there is no alternative explanation to this day). In his 1948 
paper Dirac says: "// one supposes that a particle with a single magnetic pole 
can exist and that it interacts with charged particles, the laws of quantum 
mechanics lead to the requirement that the electric charges shall be quantized 
- all charges must be integral multiples of a unit charge e connected with 
the pole strength by the formula eg = ^hc. Since electric charges are known 
to be quantized and no reason for this has yet been proposed apart from 
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the existence of magnetic poles, we have here a reason for taking magnetic 
monopoles seriously". He then goes on to say that the fact that they have 
not yet been observed may be ascribed to the large value of the quantum 
of the pole. 

- a large class of theories that include electromagnetism as a subset 
predict magnetic monopoles as solitons, as was shown by 't Hooft and 
Polyakov in 1974 [6] , and 

- if magnetic monopoles exist, Maxwell's equations are symmetric un- 
der the exchange of electric and magnetic fields. This duality symmetry 
relates small electric charge to large magnetic charge and viceversa. A gen- 
eralization of this symmetry to non-abclian theories would mean that the 
dual theory (of weakly coupled monopoles) could be used to understand 
strongly coupled non-abelian gauge theories and, in particular, confine- 
ment, for which there is no other analytic approach. While the idea is not 
new, some of the most important developments in this area are fairly re- 
cent; but they fall outside the scope of these lectures, and I refer the reader 
to an excellent review by Harvey [7]. 

In what follows 1 will take c = h = 1. 1 will also depart from Dirac's 
notation and use Qm to refer to the magnetic charge; g will be the SU (2) 
coupling constant. 

One final comment. I think everyone attending this school is aware of 
the language problems between high energy and condensed matter physi- 
cists. A relatively common source of confusion is the use of the word gauge 
symmetry, which can mean different things to the two communities. We 
all agree that electromagnetism has a gauge symmetry, it is the symmetry 
that allows local (that is, position-dependent) changes in the phase of the 
wave function and a compensating change of gauge in the vector potential 



The electromagnetic tensor F^i, = V^Ai, — V;/^^, also known as the field 
strength, is unchanged by this transformation, while the covariant derivative 
of tp, Vy^iV' — ieA^il), transforms in the same way as -0 itself (thus the 
name covariant). 

High energy physicists use the term "gauge symmetry" to indicate any 
symmetry which is local, whether or not it corresponds to a L''(l) trans- 
formation. Most condensed matter physicists, on the other hand, will talk 
about gauge symmetries to indicate a change in the phase of the wave 
function (a U{1) transformation), whether or not there are vector poten- 
tials around. Thus, the transformation 



A^{t, x) ^ A^{t, x) + V^xit. x) 



^ = 0,1,2,3 . 



(2) 



a = const 



(3) 
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is a gauge transformation in the condensed matter literature, but not in 
the high energy hterature (where it would be called a global U{1) transfor- 
mation). On the other hand a high energy physicist will talk about e.g. a 
SU{2) or 50(3) gauge transformation, meaning what is best described as 
a "position-dependent rotation" in internal space. Since there are three de- 
grees of freedom associated with rotations, we need three vector potentials 

(1) 

where r = (T^,r^,r^) are the Pauli spin matrices. The transformation law 
for the gauge potentials is 

W^(t,x) ^ M-\t,x)Wf,{t,x)Mit,x) + -M-\t,x)Vf,Mit,x) (5) 

g 

where M(t, x) is a SU (2) matrix at each point in spacetimc; note that if 
M were independent of position, as in the case of spin, there would be no 
need for gauge potentials and the symmetry would be called a global SU (2) 
symmetry. 

The transformation is non-abelian and moreover the field strength, 
which has to be generalized from F^j, = d^A,y — d^A^ to 

G^u = 5^W, - 5,W^ + 5[W^W, - W,W^] (6) 

or 

G^u = d^W_, - d,W^ + gW^ X (7) 

(with the cross product taken in internal space) is no longer invariant but 
changes like 

G^u-r^ M-\t, x) G^, ■ T Mit, x) (8) 

The transformation law of the scalars depends on the group representa- 
tion to which they belong and we list here two that will be relevant later. 
The fundamental representation of SU{2) is a doublet of complex fields 
$ = (j)2) whose transformation law and covariant derivative are 

$ ^ M~^(t, x)$ and D^^ = + gW^^ respectively. (9) 

The adjoint representation is a triplet of real scalars ^ = {(j?- ^cjP' ^(f') which, 
like the gauge potentials, can be assembled into a matrix ^ = (j) ■ t with 
transformation law 

^it,x) ^ M-\t,x)^{t,x)M{t,x) (10) 
and covariant derivative 

D^^ = V^* + 5[W^* - *W^] (11) 
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or, equivalently, 

D^^ = V^^ + gW^x^. (12) 

In these expressions, ji = 0,1,2,3. In what foUows, we wiU "work in 
temporal gauge", setting the time component of the gauge fields to zero. 
Thus, the vector potential (or gauge potential) will be a three vector and we 
will use the notation V x W + gW x W for the field strength; the expressions 
above should serve to clarify whether the cross product is taken in internal 
space, in real space or in both. 

3. Do-it-yourself magnetic monopoles 

This section is "adapted" (i.e. taken) from Coleman's 1974 Erice lectures 
[8]. He refers to this as "the monopolc hoax"', a joke to be played (or at 
least attempted) by a cunning theorist on a gullible experimenter. In view 
of the number of cunning experimenters in this audience I will refrain from 
comments and just describe here how to build your own magnetic monopole. 

1) Take a solenoid. It has to be very long and very thin so as to be 
invisible; as Coleman says, it helps if the solenoid is many miles long and 
considerably thinner than a fermi (this is very much a gedanken hoax). 

2) Put one end at the experimenter's laboratory, 

3) hide the other end, and 

4) turn on the current. 

For a gullible theorist this may pass as a magnetic monopolc, but of 
course there is a way in which the solenoid could be detected: through 
Aharonov-Bohm scattering. 

The interference pattern in a double-slit experiment is shifted when a 
solenoid is placed between the slits and the screen. Even if the particle 
trajectories remain well outside the solenoid, their wave functions ipi, ip2 
acquire a phase exp [ie f A-dl] (with e the electric charge of the particle and 
the integral taken along the particle's path); if the paths are on either side 
of the solenoid, the interference pattern changes because the probability 
amplitude + V'2p becomes 

I ie r A-dl I I ie f„ Adl i \2 i ; . ie S A-dl i |2 f-io\ 

I e Ji ipi + e V2I =\Vi + e J ip2\ , (13) 

where § A ■ dl, taken around the solenoid, measures its magnetic flux. 

Notice that, if the flux is an integer multiple of 27r/e, the solenoid be- 
comes undetectable even quantum mechanically in our gedanken experi- 
ment - it is called a Dime string. 

The vector potential of a monopole whose Dirac string is along the 
negative z-axis is given (in spherical coordinates (r, 9, (p) centred on the 
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monopole) by 



An ■ dx = y^(l — cos 9)d(p, 



(14) 



giving rise to a radial magnetic flux 



B -dS = ^ sin edOdip 



or 




47r 



(15) 



where Qm is the magnetic charge. The vector potential is singular on the 
Dirac string, which is located at ^ = tt, but regular everywhere else. 

Since the magnetic flux of the monopole is supplied by the Dirac string, 
the condition that the flux through the string should be a multiple of 
27r/e gives rise to the famous Dirac quantization condition for the mag- 
netic charge Qm of the monopole: 

= or, reintroducing *a.dc, ^ = ^(M (16) 

e 47r 2 

Thus, the existence of one monopole would be enough to force electric 
charge to be quantized! 

Note that we can use gauge invariance to change the position of the 
Dirac string; an equivalent description of this monopole is given by the 
vector potential 



You should not worry about these singularities - in the next section we will 
eliminate the Dirac string altogether. 

4. The Wu-Yang construction of Dirac monopoles 

Electromagnetism is not a theory of gauge potentials per se, but rather of 
equivalence classes of gauge potentials. This was exploited by Wu and Yang 
[9] to give a non-singular description of magnetic monopoles by "patching 
up" vector potentials that are regular in different regions, provided they 
are equivalent on the overlaps. 

Consider any sphere with non-zero radius surrounding the monopole. 
In the northern hemisphere 9 G [0,7r/2 + e] take A = Aj^ and in the 
southern hemisphere 9 G [7r/2 — e,7r] take A = As- In the overlap region. 




^Note that B does not include the singular contribution from the Dirac string. 
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7r/2 — e<^<7r/2 + e, the two descriptions are related by a regular gauge 
transformation, 

(An — As) ■ dx = Vx • dx where x{^) = tt^ ^ ■ (18) 

ZTT 

This object has magnetic charge as can be seen by computing the 
magnetic flux through the two hemispheres (using Stokes' theorem). Note 
that, even though x is not singlevalued, Vx is, and therefore the gauge 
transformation is weU defined on the gauge potentials. Moreover, single- 
valuedncss of the gauge transformation on the wave functions, e''^^, in the 
overlap equatorial region implies the Dirac quantisation condition! q^yi = 
X(27r) - x(0) = 27rAr/e. 

We have eliminated the need for a Dirac string - note that the only 
singularity in this description is at the origin, r = 0. But this must be a 
real singularity because the energy of the monopole diverges as r — >^ due 
to the 1/r^ behaviour of the magnetic field, 

E = Jd^xli\E\^ + \B\^)r. JJ"^ (19) 

Electromagnetism is not the only force in nature. There are the weak 
and strong nuclear forces, and also gravity. We think that forces may 
become unified into one kind of interaction (a Unified Theory) at high 
energies^. When theorists started to investigate possible unified theories 
they found a surprise... 



5. 't Hooft— Polyakov monopoles 

One of the first attempts to unify the electromagnetic and weak interactions 
was the 0(3) Georgi-Glashow model [10] in which the fundamental fields 
are a triplet of (real) scalars and a triplet of gauge potentials. 

^=(^^0^0^) 

The most important aspect from our point of view is the energy, since 
we are seeking to remove the divergence at r = 0. I am making many 
simplifying assumptions here (no time dependence, no electric fields), and 
only writing the terms in the energy that are relevant for the argument: 



E = j d^x [{V^ + gWx^f + CV xW + gWxWf + Xi^-^-rjY] (21) 

■^A special class are so-called Grand Unified Theories, or GUTs, where the electro- 
magnetic, weak and strong interactions are described by a single simple group. 
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First of all, the Georgi-Glashow model includes electromagnetism as a 
subset: the configuration 

^ J ! (22) 

W^ = W^ = , = A , 

where A is any solution to Maxwell's equations, is also a solution of the full 
non-abelian field equations of the Georgi-Glashow model. In particular, 
the Wu-Yang (or Dirac) monopole is a solution. However, it is not a stable 
solution! 

Indeed, the l/r^ divergence in the monopole energy is now coming from 
V X in the /(V x + gW^ x W'^Y term and could be controlled if 

and W"^ acquired non-zero values ~ 1/r. This is consistent with the fact 
that W"^ = {W^ ^ iW^) / are charged fields (the W-bosons) with charge 
zizg respectively and a magnetic moment igW~ x which couples to the 

magnetic field, so their presence can reduce the magnetic energy. On the 
other hand, such "W-condensation" has two immediate effects: one is an 
increase in energy coming from the new, non-zero (V x W^''^±gW^ x W^'^)^ 

— * 2 1 

terms; the other is that the scalar gradients ~ a W ' x (f)^ now 

diverge as 1/r. However this problem is eliminated if ~ r as r ^ 0. The 
condition (j){r = 0) = imposes a penalty in energy from the / X{(p ■ cj) — r]"^)"^ 
term, but this is finite - thus, the result is always energetically favourable 
to the singular abelian monopole that we started with. 

In the case when magnetic charge is two Dirac units, qm = ^i^/g (see 
below for an explanation of this condition), this instability leads to the 
't Hooft-Polyakov monopole [6] , a spherically symmetric configuration de- 
scribing a non-singular magnetic monopole of finite mass [11], 



(/.i = (/)^ = , (f)-^ = r]p{r) 

#1 - iW^ = iZMe^<^( _ i0\^ #3 ^ _ j.os^)_L^^ ^ 

g r smO g rsmO 

(23) 

with /(r) 1, /'(r) ~ — r and p(r) ~ r as r — > and f{r) — > 0, p(r) — >■ 1 
as r ^ oo. 

Note that only the small r behaviour of the fields has changed; in par- 
ticular, the magnetic charge of the monopole remains the same. After a 
(singular) gauge transformation it reduces to the more familiar form [6] 

(f) = p{r)r 

_(/(r-)- l)f X Vr , ^ ' 
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which shows that the 't Hooft-Polyakov monopole is a topological defect 
(usually called a hedgehog because of the way the scalar field points radially 
outwards). The zero value of the scalar field at the origin r = is forced 
by the non-trivial winding of the scalar field. 

It remains to explain why the restriction to two units of magnetic charge. 
In the Wu-Yang construction we started with a sphere of non-zero radius, 
say R, divided into two hemispheres overlapping at the equator. If the 
monopole has A'^ units of magnetic charge the gauge transformation in 
the overlap region is a phase rotation by 27riV. W-condensation replaces 
the singularity at the origin by an everywhere regular core. Since nothing 
changes outside the core, the patching condition for the 't Hooft- Polyakov 
monopole remains a 2TrN rotation for all R > rcore- But this cannot be 
true inside the core: since there is no singularity, the gauge transformation 
must also change continuously so that it becomes the identity when we 
reach r = 0. If the gauge group is U{1), this is simply not possible, and 
all monopoles are singular. But in SU{2), a 27r rotation is not continuously 
connected to the identity whereas a 4tt rotation is! Thus, only monopoles 
with even can be non-singular. Of all these, it turns out that only N = 2 
remains spherically symmetric after W-condensation. 

The existence of magnetic monopoles is a very generic prediction for a 
large class of theories containing electromagnetism. Moreover, they should 
be produced in large numbers in the early Universe [12]. The fact that we 
do not observe those monopoles is a serious challenge to cosmologists, and 
has become known as the monopole problem. 

6. Magnetic monopoles in the Weinberg— Salam model; electroweak 
strings and dumbells 

The standard model of electroweak interactions has a SU{2) x U{1) gauge 
symmetry, corresponding to weak isospin and hypcrcharge respectively. Its 
bosonic sector comprises a neutral scalar field (j)^, a charged scalar field 
0+, and the vector potentials corresponding to the massless photon A and 
three massive vector bosons: the charged W-bosons (W^) and the neutral 
Z. The fermionic sector consists of the three families of quarks and leptons 
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(25) 



In the Weinberg-Salam model, the electromagnetic and Z-fields are com- 
binations of the SU{2) and ^7(1) gauge potentials {W_ and Y respectively): 

Z = cos Oy^W^ - sin Oy^Y , = sin Oy^W^ + cos Oy^Y , (26) 



200 



where 6yu is called the weak mixing, or Weinberg, angle. Its measured value 
is sin^ Oyj K 0.23. 

In this case, the field that satisfies Maxwells equations at low energy is 
A and, being massless, it is the only vector potential that can give rise to 
long-range electric and magnetic fields and thus magnetic monopoles. Its 
configuration far from the monopole will be exactly like what was discussed 
in sections 3 and 4. Very close to the monopole, though, we would expect 
other fields to condense due to the intense magnetic field, changing the core 
structure. 

But there is a problem. Note that, since the electromagnetic field has a 
hypercharge component and hypercharge is an abelian field, isolated mag- 
netic monopoles are always singular at the origin. We are back to square 
one! 

In order to have monopoles with regular cores one has to embed the 
SU (2) X C/(l) symmetry of the Weinberg-Salam model into larger symmetry 
groups. We already mentioned Grand Unified Theories (GUTs, for short), 
where one simple group not only contains the electroweak interaction, but 
also the strong interaction. These monopoles are very heavy, because the 
unification of these forces occurs at very high energies {'^ 10^^ GeV) and 
the fields that condense at the core are very massive. Far too heavy to be 
produced in a particle accelerator. 

But there are also lighter magnetic monopoles in the Weinberg-Salam 
model: they occur as monopole-antimonopole pairs connected by a vortex 
(the vortex carries magnetic flux of the .Z^-boson, and it is usually called a 
Z-string or an electroweak string) . Such configurations were called dumbells 
by Nambu, who first considered them in 1977 [13]. The dumbell is rotating 
to avoid longitudinal collapse, and its mass is estimated at a few TeV. 

Their internal structure is rather interesting. The SU(2) fields are those 
of a 't Hooft-Polyakov monopole-antimonopole pair, while the hypercharge 
U{1) field configuration resembles that of a solenoid joining the monopole 
and antimonopole. As a result, the combination inside the solenoid is pre- 
cisely the magnetic part of the Z field, whereas the magnetic field that 
emanates from the solenoid ends is the massless electromagnetic field, and 
there are no singularities anywhere. 

In some respects, the structure of the Z-string is similar to that of a 
magnetic vortex in an Abrikosov lattice that appears in a type II super- 
conductor subjected to an external magnetic field. In the Weinberg-Salam 
model, the role of the vector potential is taken by the Z-field, and the 
order parameter is the neutral Higgs field (f)^. But there is a very impor- 
tant difference: electroweak strings are non-topological, and therefore not 
necessarily stable. 
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A review of electroweak strings where the stability issue is discussed 
in some detail can be found in [14]. Stability depends on the value of the 
Weinberg angle and the masses of the various fields. In the physical range of 
these parameters, it is found that infinitely long, straight, bare strings are 
classically unstable (they are only stable for ~ 7r/2 and mniggs < mz)- 
On the other hand, stability improves for short segments and in the presence 
of magnetic fields. Finally, the strings are superconducting and fermion 
modes on the string might also stabilize them, although this is still under 
discussion. The stability of dumbells remains an open question. 

Even in the region of parameter space where strings are (classically) sta- 
ble, the fact that stability is not topological immediately raises the question 
of whether a network of strings would form in a phase transition via the 
Kibble mechanism. We now turn to this question, and we will focus on 
the (unphysical) limit of the Weinberg -Salam model in which the SU{2) 
symmetry becomes global, while keeping the hypercharge U{1) symmetry 
local. In this limit, known as the semilocal model, electroweak strings are 
classically stable. 

7. Semilocal strings 

The semilocal model is obtained when the complex scalar field in the 
Abelian Higgs (or Landau-Ginzburg) model is replaced by an SU (2) dou- 
blet of complex fields {(f), if)) . ^ 

It is also (the bosonic sector of) the Weinberg-Salam model in the 
limit in which the SU(2) gauge coupling is set to zero; in this limit, the 
I^-bosons and the photon decouple and the symmetry is S'C/(2)giobai x 
C/(l)local- The only gauge field is the neutral Z boson, which coincides with 
the hypercharge field. The model has vortices (semilocal strings) whose 
properties are intermediate between electroweak strings and Abrikosov- 
Nielsen-Olesen vortices (see [14]). 

The energy per unit length of cylindrically symmetric configurations is 



and it turns out that, even though the vacuum manifold is a three-sphere. 



which is simply connected, 7ri(S'^) = 1, there are stable strings if the scalar 
mass is smaller than the vector mass. These strings are not only stable 



E = 




1.2 , i2 ' 

^ +"0 = Y 



(28) 



^Similar systems have been considered in the condensed matter literature [15, 16]. 



202 



classically, they are also stable to semiclassical tunnelling and to breaking 
by monopole pairs. If mgcaiar > JTivector the strings are classically unstable, 
and if the masses arc equal there is a two-parameter family of configurations 
with the same energy where the quantised magnetic flux 27r/gf spreads over 
an arbitrarily large core width. 

Semilocal strings, like electroweak strings, can have open ends; but the 
monopoles at the ends are global monopoles and have a long-range inter- 
action. 

The stability of semilocal strings is well understood in terms of the 
competition between gradient energy and potential energy. The gauge field 

can compensate gradients in the (complex) phase of either of the two scalar 
fields, but it can only compensate both gradients simultaneously if there is a 
correlation between the phases of the two scalars; this correlation is present 
in the string solution, which optimizes gradient energy, but at the expense 
of a large concentration of potential energy at the core. If X/q'^ is very large, 
it becomes energetically favourable to break the phase correlation in order 
to reduce potential energy, and the string is destroyed. If A is small, the cost 
in gradient energy to dissolve the string becomes too high and the string is 
stable. 

Since the strings are non-topological, the question immediately arises 
as to whether they would form at all in a phase transition. In particular, 
stability depends on the mass ratio m^scalar /'^vector SLud we would expect 
the formation rates to reflect this dependence. 

We now turn to a first attempt to answer these questions through nu- 
merical simulations. Note that setting ijj = in the semilocal model gives 
the Abelian Higgs (or Landau-Ginzburg) model, thus making it possible to 
compare the rates of formation of semilocal strings to those of topological 
strings in a similar environment. The conclusion seems to be that semilocal 
strings do form, and in some cases with number densities comparable to 
those of their topological counterparts. 

8. Numerical simulations of semilocal string formation 

For details of the simulations we refer the reader to [17, 18]. Here we will 
just point out the main features and results, summarized in figure 5. 

Space is discretized into a lattice with periodic boundary conditions. 
The equations of motion are solved numerically using a standard staggered 
leapfrog method, and a dissipation term (rjcp, rj^ or rjYi) is added to each 
equation to reduce the relaxation time. A range of strengths of dissipation 
was tested, and it did not significantly affect the number densities obtained; 
the simulations displayed in the figures all have have rj = 0.5. Note that in 
an expanding Universe the expansion rate would act as a sort of viscosity. 
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Figure 1. The flux tube structure in a two-dimensional semilocal string simulation with 
/3 = 0.05. The upper panel {t = 0) shows the initial condition after the process described 
in the text. The lower panel shows the configuration resolved into five flux tubes by a 
short period of dynamical evolution [t = 100). These flux tubes are semilocal vortices. 
Note the diflferent numbers of upward and downward pointing flux tubes, despite the 
zero net flux boundary condition. The missing flux resides in the smaller 'nodules', made 
long-lived by the numerical viscosity; the expansion of the universe could have a similar 
effect and preserve these 'skyrmionic' configurations [20]. 



though 77 would typically not be constant. However, this is not meant to be 
a cosmological simulation since spacetime is flat. 

We work in temporal gauge. Then Gauss' law becomes a constraint 
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derived from the gauge choice Yq = 0, and is used to test the stabihty of 
the code. 

The inverse vector mass is taken as the unit of length and time, and rj^ 
as the unit of energy. In these units the dynamics is governed by a single 
parameter, 

2A 

/3 = "^-scalar/ "^vector = (29) 

which also determines the stability of the straight string solutions: they are 
stable if /3 < 1 and unstable if /3 > 1. ^ 

The number density of defects is estimated by a multi-step process. 
The initial conditions are obtained with a generalization to non-topological 
strings of what is known to cosmologists as the Vachaspati-Vilenkin algo- 
rithm [19], followed by a short period of dynamical evolution. 

• First, we generate a random initial configuration for the scalar fields 
drawn from the vacuum manifold, which is not discretised. If space is a grid 
of dimension iV^, the correlation length is chosen to be some number p of 
grid points {p = 16 in [17, 18]; the size of the lattice is either = 64 or 
N = 256). To obtain a reasonably smooth configuration for the scalar fields, 
we assign random vacuum values on a (N/p)^ subgrid and interpolate the 
scalar field smoothly onto the full grid. 

• We then find the gauge field configuration that minimizes the energy 
in this fixed scalar background. 

Two dimensional test simulations have shown that the energy minimiza- 
tion is redundant, since the early stages of dynamical evolution carry out 
this role anyway; for simplicity, we used in practice a gauge field configu- 
ration which was close, but not equal, to the real minimum. 

• An example of the initial conditions generated with this algorithm 
can be seen in Figure 1 in the case of a two-dimensional toy model with 
translational invariance in one dimension, say z. The plot shows magnetic 
field on the x — y plane, which has a complicated flux structure with extrema 
of different values (top panel of Fig. 1), and it is far from clear which of 
these, if any, might evolve to form semilocal vortices; in order to resolve 
this ambiguity, the initial configurations are evolved forward in time with 
zero initial velocities for the fields. After a short transient, in the unstable 
regime /3 > 1 the fiux quickly dissipates leaving no strings. By contrast, 
in the stable regime /3 < 1 stringlike features emerge when configurations 
in the "basin of attraction" of the semilocal string relax unambiguously 
into vortices (bottom panel of Fig. 1). We only count vortices after this 
relaxation process. 

is also the parameter that distinguishes type I (/? < 1) from type II (/? > 1) 
superconductors. 
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Figure 2. Loop formation from semilocal string segments. The figure shows two snap- 
shots, at f = 70 and t = 80, of a 64'^ numerical simulation of a network of semilocal 
strings with /3 = 0.05, where the ends of an open segment of string join up to form a 
closed loop. Subsequently the loops seem to behave like those of topological cosmic string, 
contracting and disappearing. 



• One important point is that strings are always identified with the 
location of magnetic flux tubes, rather than by the zeroes of the scalar 
field. Figure 3 shows two snapshots of a simulation on a 256^ lattice. In 
order to make the figure, and also to compute the number density, we need 
to set a magnetic flux threshold - the strings that can be seen in that 
figure are made of those points in which the magnetic field exceeds half 
the maximum value in an Abrikosov-Nielsen-Olesen vortex. After the short 
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transient, this threshold can be modified without significantly affecting the 
results, which are ratios of semilocal to topological string number densi- 
ties (shown in figure 5). The error bars include, among other things, this 
threshold dependence as well as the dispersion between runs. 

Since the initial conditions are somewhat artificial, the results were 
checked against various other choices of initial conditions, in particular 
different initial conditions for the gauge field and also initial conditions for 
the scalar field closer to a thermal environment (see Fig. 4). However, it 
cannot be sufficiently stressed that it is not a realistic thermal simulation of 
the phase transition. There is no thermal noise and all the initial conditions 
in [17, 18] had zero initial velocities for the fields - initial conditions with 
non-zero field momenta have not yet been investigated -. What the simula- 
tions show is that the rate of formation of these non-topological vortices is 
not only non-zero but in certain cases it could even be comparable to that 
of their topological counterparts. Obviously the details of the transition can 
and should now be investigated in full. 

9. Discussion and outlook 

The physical mechanism behind the observed string formation by accretion 

of magnetic field in Figure 1 or by growth of string segments in Figures 
2 and 3 is not very different from the formation of Abrikosov lattices in a 
type II superconductor, even though in a system such as the early Universe 
it makes no sense to talk about external magnetic fields. Because of the 
y^<I>^ term in the energy, the magnetic field does not like to coexist with 
the superconducting (0 / 0) phase of the scalars. 

In the topological strings case one usually argues that the non-zero 
winding in the phase of the scalar field forces a zero of the Higgs, and 
magnetic flux gathers there in a vortex. Conversely, what we are observing 
in the semilocal model is that if there is a sufficiently large concentration 
of magnetic flux in a small region, for instance near the end of a string 
segment or maybe due to a fluctuation, a line of zeroes of the Higgs can 
develop there. Small segments of string, on the other hand, will tend to 
contract and disappear. 

It seems clear that these results should extend to electroweak strings 
with g ^ as long as we remain in the region of parameters where the 
strings are classically stable ^. In the region of stability, we expect a non-zero 
density of electroweak vortices to form in a phase transition. Preliminary 

®Note that the monopoles at the ends of semilocal strings are global, whereas those at 
g ^ will have finite size cores. However, the core size grows as 1/g and, for sufficiently 
small g, it may be larger than the average distance between string segments, causing the 
strings to grow 
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Figure 3. The growth of string segments to form longer strings. The figure shows two 
snapshots, at time t = QQ and t = 70 of a large 256^ numerical simulation of a network 
of semilocal strings with /3 = 0.05. Note several joinings of string segments, e.g. two 
separate joinings on the long central string, and the disappearance of some loops. The 
different apparent thickness of strings is entirely an effect of perspective. The simulation 
was performed on the Cray T3E at the National Energy Research Scientific Computing 
Center (NERSC) in Berkeley. 

results seem to confirm this picture [21], so obviously this problem deserves 
attention. 

We may be a long way from understanding the formation of magnetic 
monopoles in a phase transition, but it is possible that particle accelerators 
will show signatures of monopole-antimonopole pairs in the not too distant 
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Figure 4- The number of somilocal strings formed per initial two-dimensional correlation 
volume in a toy model with translational invarianco in one direction. Each point is an av- 
erage over ten sinmlations. Squares indicate that the vacuum initial conditions described 
in the text were used, while open circles indicate that non-vacuum (more thermal) initial 
conditions were used. 




Figure 5. The ratio of lengths of semilocal and cosmic strings. 



future. If Nambu's prediction is correct, dumbclls could be the first soHton- 
like objects in the standard model of particle physics to be observed. 
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